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ii to solve the non-linear cquations of motion of a generalized VTOL 

ft aircraft taking off from or landing on the deck of a ship moveny, in an 
i irregular or random seaway. 
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This version of LARC-I is limited to longitudinal nsttons 
but is designed for eventual expansion to all degices of freedom. 

The LARC-I programs makes use of ship motion amplitudes and 
: frequencies derived separately in a ship motions program, wherein 
the forcing functions of the seaway are based on a stochastic representa}- 
tion of the waves for any given sea state. The pitching and heaving | 
motions of the ship are transmitted to the aircraft by a realistic 
simulation of the landing gear. . 

The program permits solution of an arbitrary maneuver, defined 
by a sequence of aircraft contxol changes, and composed of any number 
oF segments of arbitrary length, in each of which the controls remain 

fined. . 
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fhe equations of mution for a gencralized VTO!, aircraft are 
wertleularized by adding a force and moment module for each specific 
aLrerart. This version of LARC-I contains one such module for the 
“E-S3A, D helicopters. 

Dynamic validation with flight test data is considered excellent 
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SUMMARY 
Project Description 
This report describes the work accomplished by CADCOM, 


Inc., under O“R Contract N00014-72-C-0531, titled “Sea Control 


Ship/sijcraft Motions". The work was sponsored jointly by 


une Offien of Naval Research and the Naval Air Systems 
Command, 

The purpos:. of the project was to develop an analytical, 
computer~ augmented, simulation to evnciaee: (a) a launch/ 


recovery effectiveness evalvation of ship-aireraft combina- 


tions under a variety of wind/wave and ship motion conditions, 


and (b) a specific effectiveness evaluation directed to a 


j 

Sea Control Ship and the VTOL aircraft to be associated with : 
its operation. The principal objective was the construction ey ; 
and subsequent avaluatian of an accurate and representative vv : 
simulation model, : 
Tt was established that certain limiting aasunptions 3 

and sinplifications would be necessary in order to make 3 
the development of a simulation posgsibie under the given 
funding and time constraints and that the exercise of the sedel 
would be necessarily confined to two candidate sea control : 
ships and seo aireraft (Cd-53 and Harrier). - : 


the eimslation package which was developed under the ‘ 
contract ia called “LANC-I" {LAunch ard Recovery Compatibility}. 
LARC=! situlates the motion of an aireratt in an irregular 
Seaway. ALthaugh the motions of | ¢ aiteratt are con- 


fined ty the x-z plane and those of the shig to heave, 


rey 
is 4 


pitch and a forward motion, the LARC-I program is 


ax 


Psp 


structured to allow an orderly modular expansion into 
a version which will sinuiate the motion in all six degrees 
of freedom, longitudinal and lateral. | an 

The ship motions selected for the exercise o. the model 
are parametric in nature and hence can represent the motions 
‘f a number of candidate air-capable ships. The aireraft 
used in the model exercine and generation of data wa: che 
CH-53D helicopter. During the performance period, 4DCOM, 
ONR, and NAVAIR attempted to obtain data on the |. .vier 
aircraft without success. Thus, the latter ai craft vas 
eliminated from the computational phage of th: project and, 


instead, additional data was generated for the CH-53D, . . 
Conclusions 


Data on the Harriex aircraft was received by CADCOM 
after the completion af this project, and it is included in 
this repert (Appendix £) for reference purposes, 

AS a reault ef ti.is wark, the tentative conclusion 
haa been reached thst pitching and heaving platfors motions 
wiki not limit normal take-off and resovery maneuvers of 2 
VTOL @irerait of this type. Natural frequeacias ars too 
far apart to couple. Although thera ia an inersase in 
control frequencies due tao fuselege pitching momanta, the 
rejuiread pilot rewponse ie within normal limits, Sreeifically: 

&. For the ship motian tomlitiona investigated, 16 
serious stability problem ocsurred which would Lisit the 


take-otf capability of the CH-$1) nalicopter. The static 
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/ 
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and dynamic stability margins of the aircraft are sufficient 
to allow safe take-off maneuvers with acceptable levels of 


pilot effort. 


b. In no caSe were landing gear loads in excess of 


the design limit loads as given in References 6 and 7. (Note) 


Recommendations 


x Asa result of its experience with surveying the state 
a of the art, developing the mathematical model, constructing 
the computer simulation, and applying LARC-I to realistic 
ship/aircraft operations. CADCOM makes the Following recom- 
iq inendations: 

4 : a. Successful and realistic general use of the LARC 

aa simulation as an analytical tool in system analysis will 
4 require the extension of the progrem to all six degrees of . 
q freedom for both ship and aircraft. It is apparent from 
the results of the longitudinal motien studies that high~ 
frequency cyclic response will result in roll farces which 
may couple with landing gear forces in the roll plane. 

b. Real-time or quasgi-real-time flight control should 


be added to the simulation to allow the operator to “fly* the 


: ; = 2 
Oa . ae =: ea 

See eUey ATE. Coa PLAC R Te 
Sls So as OD hes ii Bee A Sa RYO 


e: model. His reference system would be provided by an inter- 


sk 


4 active couputer graphica display. In the present version, 
the model operates in short (5-8 second) control-fixed seg~ 
ments. In the proposed interactively controlled model, the 
Notes oasies linitt loads « main teers syaaa? 


(each gear - cond. *. 5, LL) 


4 nose gear: 173008 
4 {each gear ~ cond. L. L. 3 PT.) 
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high sampling rate control loop would allow adaptive control 


in a real-time mode which sore accurately reflects pilot 


response, | 

ec. Investigation of the mathematical representation of e 
the random sea with the intention of resolving the problems | 
discussed in Appendix D should be continued. 
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I. INTRODUCTION 

Frequently the effect of the motion of a ship in 
moderate to heavy seas on the aircraft and/or the ship 
design and performance characteristics necessary for a 
high probability for a successful launch or recovery oz 
aircraft have been given little consideration during the 
design phase. iin fact, the design trade-offs between 
the ship and the aircraft performance characteristics are 
probably unknown to a large degree. 

Recently the Sea Control Ship has emerged as a signifi- 
cant advancement ir our defense posture, As presenily con-~ 
ceived, the Sea Control Ship is a small, versatile combatant 
with significant air capability. If this new concept is to 


succeed, significant effort must be expended *o insure that 


en 


small ships can serve as adequate olatferms for aircraft. 
ince the platforms on such ships will of necessity be 
smaller, the motions of those platforms will consequently 
bo more sevare than those of praesent aircraft carriers. 
GQou Of the first steps in the process of evaluating 
the alr capable potential of candidate ship types and/or 
Gesiqna is ta determine their seakeaping characteristics, 
Qnee the response cf the candidate ship to various sea 
atates is known, this informatian can be assessed together 
With the maneuverability and control] characteristics of 
candidate aircraft to determine the petential as a Sez 


Centrol Ship or aircraft. 
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Using Gompueersaided analysis and interactive graphics, 
techniques it is possible to simulate the dndiyiadalnoetone 
of both the aircraft and the ship and hence their relative 
motions. This leads to the ability to estimate the proba-~ 
bility of a successful launch/recovery with existing aircraft 
from existing or proposed snips. In addition, it is possibie 
to simulate the motion of proposed ship concepts and thus 
to determine the probability of a successful launch/recovery 
with existing or proposed aircraft. 

In the following discussion, the motion of the ship is 
first discussed in Section II, then the motion of the 
aircraft in Section III, and finally the combined motion as 
used in the simulation as well as the simulation itself in 


Section IV. 
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II. SHIP MOTIONS 


The characteristics of the ship which were used in 
this simulation are shown in Table II-1. Further details 
are given in detail in References 1 and 2. The determina- 
tion of the ship resronse amplitude operator (RAO) was by 
means cf the linear superposition theory described in 
References 1 through 3 except for the pitch RAO at a 
Froude number of 0.30. These responses were derived from | 
model tests. The theoretical and experimental values are 


shown in Figure II-1l. 


TABLE II-1 
MODEL M 
PARAMETER SHIP MODEL VALUES 
i2 
LENGTH BETWEEN PERPENDICULARS 
LSP, FT 12.69 
DRAFT AT OL; FT 0.725 
BLOCK COSFFICIENT 0.632 
BEAM, FT 1,93 
CENTEK OF BUOYANCY, FT. AFT OF 7% 0.208 
STATIONS 10.33 
CENTER OF “LOTATION, PY AFT OF 3B 324 
STATIONS 10,83 
RADIUS OF GYRATION/LEP, FT 0.256 


The full scale length of the ship was assumed to be 630! 
LWL and LBP. The maximum sceed of the vessel was assumad to 


be about 25 knots, corresponding to a Froude number 
vu 
Fr = « 0.30 tIi-k 
Vet 
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The data shown in Figure II-l are for Froude numbers of 


0.10, 0.20 and 0.36, and for wave lengths ranging from 


L/A4 = 0.50 to 1.10. In only one case were the experimental 


results substituted for the theoretical curve; the pitch 
response at Fr = 0.30. Except for this case, theory 
and experiment are in close agreement. At Fr = 0.30, 
the theory overestimates the pitch response by up to 3 
percent. The authors of Ref. 1 attribute this to large 
“steady state” waves which may be expected at and above 
the design hull speed. The design speed for the Model M 
hull is about Fr ~ 06.28, or for a 63G6' ship, about 23.5 
knots. Unless the interaction betwcen the oscillatory 
motion and "steady state" Kelvin wave pattern is accounted 
for in the tueory, strip theory may be expected to over~ 
estimate amplitudes in the speed regime where such waves 
affect the pitching response, For pitching amplitudes 
at Fr ow 0.30 the dashed curve of Figure IYI-le was used. 
This curve was coincicent with the theoretieal curve 
between L/A © 0,9 ta t.2. 
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The spectral valves for significant wave heights of 15, 


20, 25 and 30 f£t are shown in Figure II-2. This formula 
gives an approximation to real sea condirions which has 
been recommended by the 11th International Tcwing Tank 
Conference (1966) for use “when information on typical sea 
spectra is not available." The shapes of spectral curves 
derived from data caliected in different sea locations 

and at different periods of the vear show wide variations. 
In future studies, it is recommended that more representa- 
tive data be used. At this stage fin this project) the 
Pierson-Moskowitz spectrum .s adequate since the princi- 
pal objective of the investigations described herein is 
the demonstration of the simulation model. 

The curves of Figure II-2 approximate the spectra 
which would be derived from data collected at a stationary 
point. If the point were moving at a steady velocity in a 
etraight line, the curves would be displaced to the right 
as shown in Figure It-3. The displaced curve in Figure Ir-3 
is the 20-foot spectrum as it would appear to a ship moving 


at 25.3 knots or 42.8 FES. To map the spectral data from 


(the stationary & - y plane (where &w is wave frequersy and 


¥ i@ wave headine!:.to tae moving wT e_Plane, where we is 


the frequency of encounter, where Bey and » are related by 


cos x} Ire} 
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Figure Il-2 Tho Plorson-Hoskowits Sea Spectrum for 
AS-, 20<, 25-, and 30-Feet Significant 
have Heights 
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Spectrum for “30' Wave" with 
‘ Respect to Stationary Point 
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Figure ITl-3 Sea Spectra for StaLlioenary and 
Moving Point - Significant Nave 
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it is necessary to make use of the Jacobean of the 


transformation 
jau aw 
gtx Balu gj Se Me 11-4 
Wet Xe 3 WetXe ax ax 
IWe, 8X 


Por the head-on condition, where the ship is overtaking 

the waves and traveling in the same direction, the Jacobean 
is 

II-5 


1 
(1. + Aue y) 172 


Each ordinate of the "stationary" curve is multiplied by 


the above expression to give the corresponding ordinate 


a of the transformed spectra, as it woulé appear to a moving 
ship (Reference 3). 
ws Typical values of wave frequency, encounter frequencies 


and wave lengths for the €30' ship at Fr = 0.30 are shown 


in Table Il-z2 below. 


a PARLE 1&2 

q L/w a He we 

: 4 1575 +358 265 

q . 78? 507 356 

: 1.3 525 622 »405 

5 The wave iength \ and the wave frequency ws are related 


= by the equation 


q ; z is ? 


te ce Bea | 
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2 1/2 
We §.118 () Vr Ii-6 


For U = 42.3 fps, w ia given by 


~« 3 f3 = 
wo Fs + 1/2 tet 3a, Ii-7 


The transformed sea spectra and the ships response 
amplitude operator squared are now multipied ordinate-by- 
ordinate to obtain the response amplitude spectra for a 
given sea state and ship velocity. Examples are given 
in Figure II-4 for significant Wave heights cf 15 ~ 20 ft 
(sea state 6) and 25 ~ 39 ft (sea state 7), for the ship 
Model M moving at Fr = 0,30. Thus, the respense enercy 


spectrum aa a function of encounter fre mauence w, is 


, er eee . ae 
5, fa.) = [Ra }] 80) 1I-8 
The tatal energy of the response spectrum ia oy definition 
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and assuming a Rayleigh c'atribution of che maxima (“ampli~ 


tudes") the Zellowiag response values may be expres#ed: 
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TABLE II-3 


SIGNIFICANT MOTIONS 
(Avg. o£ 1/3 Highest) 


SIGNIFICANT HEAVE PITCH 

WAVE HEIGHT AMPLITUDE AMPLITUDE 
(FT.) (FT. ) (DEGREES) 
15 8.932 2.113 
20 12.207 3.163 
25 14,514 4.322 
30 18.035 5,492 


The amplitudes shown in Table II-3 are those used 
in the pregram to calculate motions at the landing locatic. 
The Piarson~Moskowitz Spectrum which was used in 
this determination of herve and pitch asplitudes expresses 
the energy in fully developed seas and it depands cn 
only one parameter, the significant wave height wre 
With moro parameters ane could represent the sea gnectra 
actually anceuntered which are for the most part net | 
fully developed. Or, mere accurately, a spectrum could be sgt 
used, derived from actual wave data at a given Location. 
Thit would provide the sighest measure of accuracy but 
wand necessarily apply only to sonmditions which may ovucur 
at that geographic ieeatian a8 that cise of ycoar. The 
Piergon-Noekowits apecttrum is used in thie atudy because 


it is simpic io apply ard gives reasonable predietions, 
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. : f sa 
Ordinarily, the simulation would determine the heave - 
amplitude of any point other than the center of gravity 
by evaluation of the heave response spectrum 
We eee 
(a) ] 


TR, yg MW) l? = ER, og (wd]? +L? UR 


ZX 


Fs 24R_ 


8,¢c9g 


g (wR W) Cos (8 (w,) - e(w.)) II-1ll 


7c 6,cg 


= where R, , is the response amplitude operator in heave (z) 
e 


at any point x 


R is the RAO in heave at the c.g. 
Z,Cg 


R is the RAV in pitch at the e.g. 


8,cg 


. S is the distance from the c.g. to the point x, 

®, is the encaunter frequency and 

s 6,€ are the phase angles for pitch and heave respectively, } 
both measuring the lead of the ship response ba 


with respect to the maximum wave elevation at the 
a midship location. 


The ship vertical displacement spectrum at any point other 
than the c.g. due to the combined motions of heave and pitch 
are deotarmined from this expresaion in the same fashion 


ag the heave and pitch separately at the c.g. That is: 


a? a ‘ Yon z 
Rey UR, yfogll S,(a,) II-}2 


and the significant response ammiitude is 
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Se is the RAO in the z direction (vertical) 
-"2,X at the point x due to pitch and heave at 


the c.g. 


Z,.(w ) is the significant vertical displacement 
at the point x. 


Although the results derived from this procedure are statis- 
tically correct, they are difficult to use because 

the motions due to pitch and heave are coupled in the 
resulting amplitude Z(w,). Use of this amplitude alone 
precludes the analysis of the effects of heave and pitch 
motions at the point x separately. It would be preferable 


to be able to use an expression of the form 


( : , ‘ b—3 ci - : - 
ig 3, (w.) = 2, cos (wit - 6) + £6,008 (wt - €) TI-14 


g 


whese Zz is the vertical displacement at x 
Seg is the maximum vertical displacementto the og 


ee is the maximum pltch displacement, 
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Oc more convaniantly 
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where ve «4 fot = §} 1X-15 
and 6° « (44 « } 
At wt # 0 cos (wt) “Jf 
sin(w t) x 0 


and 8) = 3. + {0 egSin 8” 1i-26 


eg 
If-12 


Sera ot - 
Rit deca Gaede a ee et ea mish bcatwe Rea Rei Aches SE ek cat Enh: 


Cee a lee Wak Stach aes, AERA Aen es a 


Rati. : ae Stars 
A as Scene eet Let OA 


If this is set equal to % (0), the coupled pitch/heave 
deflection at x, cos 6° and hence 6‘ can be determined. 
An example of a calculation is shown below to illus- 
trate the use of this method. 
The ship motions program predicts the following 
motions at stations 10 (approximate c.g.) and 5 (157° 


forward) for a 30' wave condition 
nem STATION os 


10 5 
VERTICAL 18.03° 28.56' 
PITCH 5.49° 5.49° 


(6°) #6’ = 28:53. - 18.03 . 9,698 Rap. 2-17 
157 x 0.0958 


At first this result appears to contradict the phase re-~ 
lations predicted by the ship motions program. Figure [I-5 
shows these relations for the L/A range from 0.4 to 1.2 

for the Model M. The heave and pitch motions are everywhere 
out of phase by between 90° and 135°. Rowever since the 
pitch always laga the heave by this amount, the vertical 
motions are additive and thus Equation II-15 is a very good 
approximation of the inaximum absolute deflection when & 

in determined in this manner. The harmonic motien at tho 
point x can be assume] very noarly equal to 


2 =f! cos {wt} 
ss - bs yr-18 


‘.* Beg coe fof) 
with 2° determined by Eg. 1i-16 and with 6 caiculated 


&s in Eq. [1-17 fer the particular point 157" ahead of the 


I-33 


Ee W Tepoy,, toy sorbuy 
: Sseud Te AusUTAedxg_ pur Teopyezosyy gry eanbht a, 


Y/T t/t t/'1 
oT 8° 9° O°T g* g° ori $° 9° 


e | i aes T ; | a meee pe eee ee ete 5 oF i EE. eg eee aN gn ee TTE yee 
cs | Sc ; } 


UE Doe te eA 


oe | 


- 
a 


“ge 


sTuttat tke baste AY, 


s, 
i 
. 
¥ 
ia 2 ARR SS cre ea ADLD oe, 


i 
sh 
FIONV TSVHE 


O 
© 
A 
apes 


4 
oe \ aera : 
: c ‘ oO’. | ; | 
= | aAVZH ao SAvan | 
\ “AD. ar ase ie oe : AD ery Ps 6 ~ : + « A 4 9 ; 
A \ | oe 0 a a 
E SAWSH y is ‘ 
: Shee hola 2 oO | tice ee ee bina ee bs a, Aapbsttese oneal enol 
Fe €°O ="°us 2°99 = "UZ i*o =~ ‘yg 


aye WR tee 
Sa sate abe Neca We wade Ka 2: 


ro 


-, 


cg, for a sea state spectrum based on a. significant 
wave height hi 73 = 30 and a ship forward velocity of 25 
knots, . 

The final problem, which must be solved in order to 
nwke the ship motions data useful in an investigation such 
as this, is the selection of a realistic value for Wes 
The procedure which has been used is to set the value of 
the encounter frequency equal to the ratio of the integrals 
ofthe amplitude and velocity response spectra in heave. 
(The difference in this value and the corresponding value 
for the pitch spectra is always very small, never greater 
than 0.02 rad/sec,) In effect, this is the ratio of the 
total energy of the response spectrum and the first moment i} 
of that spectrun, or 


a 
i 


w = I-19 


Pig rs, (w,) da, 
: da, 


e 
(It makes no difference, incidentally, whether the sea 
spectrum 8. and the RAG are computed cor wave fraquencies, 
rs or ancounter frequencies, te That is to say the 

same valuas result when the sea spectrum S, (wi and RAO, 
R(x}, arte transformed separately into the encounter fre~- 
quency plane, and then multiplied, (Rtg)? 5, (ug) to give 
the Ship response apeetram Sp (weds as when the response 
epectrus Si(s) transforms into the same plane without 


prioe shifting of the component curves.) 
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Eq. II-19 gives the value of a weighted average 
where the weights are the squared amplitudes of the ship 
response associat2d with that particular frequency. Since 
the Sp curve does not by nature possess multiple maxima, . 
Eq. II-19 computes a weighted value associated with the 
wave components (wave lengths) to which the ship can 
respond and to wave heights in the sea spectrum, §,(,)- 
which possess significant energy content to excite the ship. 
Since wave length and wave height are independent, this 


is the only feasible way of associating these values. 


cast, = 1 DPR Aaa OP LTS ITE ES RSE CARE ae SE OR OPE eee Eee 


It 4s important to note tnat We is not necessarily the 


most probable value. For example, if the ship response 


Ek RABE TO tage 


. curve is flat through a range of frequencies, thus 


aS, (w_) 
wae a =z 0 we < We < We 
dv : 2 
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i then the value computed by Eq. TI-19 will be between 

Ye, and Yeni not necessarily midway between. In this 
ease, very small changes in the amplitudes of S, (uw) 

and Rla,), while not affecting the value of 


E. fSgla)dw 


appreciably, are sufficient to effect the value of 


we to the extent that it may have any value between 
gig 

w® and. . Since the total value of E_ has not been 

e @, r 

signficantly changed, all value# of & in the range 


$1g 
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Yer te Yes are nearly equally probable. Fortunately, 
the average ship acts like a narrow band filter and 
responds only to a narrow range of wave lengths approximately 
the same length as the ship (i.e. LWL + 0.25 LWL). Since 
the period of most ship motions of interest in this study 
(say 10 - 15 seconds) is long in comparisun to most 
motions of interest in the aircraft (1 - 5 seconds) the 
value of Ye is not critical. This observation is further 
supported since Eq. II-19 assures us that the value comouted 
lies in the range of nearly equal probability. 

The values used in this study for a 3)' wave con- 


dition, 630! ship moving at 25 knots into head seas were: 


4 


heave amplitude 18.035 


12.320 {3 


n 


velocity amplitude 


w 
Pog 0.679 


it 


203 /u* 


« 


er since wave length A 


ere » 44y' = 0.70 x LWL 


Additianal results are shown in Table II-4, 
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PBL Tine SHIP CASE. 
I TI 
V+ SHIP VELOCITY, KNOTS 25 25 
Na, MAX. HEAVE ampLitupe?), pp 18.035 11.207 
Nes MAX. PLTCH ampuirupe'™?), peg 55.492 3.162 
_lgs NOH DIMENSIONAL ENCOUNTER 
reg. (2) 3.45 3.45 
5°, PHASE ANGLE OF PITCH WITH 
RESPECT TO HEAVE, RAD .698 1.357 
, : : 
i_! COORDINATES OF SHIP 
es FT 0 0 
2. 1 CG AT t = 0 
i. } 
oS: 
xX. r JAMS ( 
ip) COORDINATES OF LANDING |. (3) (3) 
LR, 
hy ,3 SIGNIFICANT WAVE HEIGHT, FT 30 20 
SS SEA STATE 7 6 


(1) Significant (i.e., average of one-third highest) 
values. In the program, the parameter is the ratio 
of the “amplitude® ny te 1/2 the “height® his: 

thus 18.053/i5=1.21. 


(2) uy tig Yhe’g where L is the ship length, 


(3) Thres locations were studied; midships and + 100° 
-from midships, all 26.5' above W.L, 
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III. AIRCRAFT MOTICNS 

As is well known, the equations of motion for an air- 
craft can be written in a variety of forms, each form more 
suitable te the solution of a specific problem (such as 
stability, response to gusts, contro), response, or automatic 
contrcl design), and each form subject to advantages and 
limitations, In the past, most stability and control work 
has been accomplished with the linearized equations of 


small disturbance theory. (Reference 4, section 4-14) , 


Since the stability of the aircraft being studied 
has already been demonstrated, Stability for normal flight 
conditions is not a primary consideration.Thus, the force 
and moment equations themselves are used in this study. 

The principal advantages in using the force and moment : {3 
equations lies in the ability to handle (a) large angular 
displacements (e.g. greater than 10°), (b) non-linear ground 

and air reactions, and (ec) the coupling reactions hetween 

longitudinal and lateral forces, 

Finally, the computer has removed most of the practical 
difficulties in the selution of the governing non-linear 
coupled equations. 

The generalized equations of motion of an aircraft, 
which aro used in this simulation, are given in Appendix 5. 
These equations are identical to the equations presented 
in Referance 4, Section 4, except for the addition of the 


force and moment cquation for the landing gear. They are 
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generalized in the sense that they apply to a configuration 
consisting of any combination of the following components: 

(a) Single rotor 

(b) Two rotors in tandem rotor configuration 

(c) Fuselage 

(d) Horizontal tail 

(e) Vertical tail 

(f) Tail rotor 

(g) Propellers or jet engine 

(h} Lift engine or deflectable thrust 

(i) Wings 


(j) Various stabilization devices 
For the present simulation, which is limited to the 
pitch plane, the governing equations of motion are: 


The X-Furce Equation: 
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(2). = Dy sin(a-e.) + Iycos (a-€,,) 
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the Pitehing Moment Equation: 
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SIGN CONVENTION 


POSITIVE 
x FWD 
Zz DOWN 
8 NOSE UP 
x FWD 
z DOWN 
8 NOSE UP 
a pes. nose up (cw rotation From wind 
vector to x-axis) 
B pos. if L.E. swash plate is down 
Y | uP 
Ci Eng fp pos, if increases local a 
Ly — FWD of CG. (in pos x-direction) 
i. BELGW CG. (in pes a-direction} 


Forecs aad momanrks ag shown in Figure TfI-k. 
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where the nctativr is given in Pigure III-1 and Apvendix 

B. The sign conventions are given in Figure III-1] and 
Table III-1. Assuming that the external forces anc moments : a: P 
are known, these eyuations represent three equations in the | 
three unknowns u, w, and 6 and can thus, in principle, be 

solved. In the present simulation this is accomplished 

mumerically using a fourtheorder Runge Kutta integration 


scheme. 

The method of determining the forces and moments of 
the governing equations for the particular aireraft considered 
in this simulation, the CH-53D, is given below. Further 
information applicable to the CH~53D aircraft is given in 


Appendix C. These methods are considered in the following 


order: 
Landing Gear Forces 
Rotox Forces, Anglo of Attack, and Engine Power 
Petar Hub Moment | 
Fusaluge Farces and Momenis 


Tail Sorces and Moments 
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intusference Angles 


Automatic Flight Control Syatem (AFCS} 
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2hG Ground Effect 
iii-i Landing Gaar Forcea 
{3} Poevsatic shock strut Isothermal compression is 


agaumed, Maxisum stroke on both geare is 12°. In the 
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following formuls, atroke ig measured from tne fully 
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compress2* position. The stroke vs. pressure charac- oo 


Seat 


teristics. (Ref. 7-9) for the landing gears are as 


follows: 
. PRESSURE ,psi 
STROKF . ' MAIN NOSE 
STATIC 2" 850 
67° EXT e" 365 
EXTENDED 12" 187 275 


Nose Gear Determine cylinder length: 


Fe (L + 4) 
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275 (L + 4) 
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12.22" from 67% EXT. 


al 
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Pressure, fully compressed 
Pe L. - Pe Le 
Po w Pe le 

Le 
Lead at any atroke 


= 273 X 16-22 . i0¢0 psi 


DD = PolpA : 
LOAD - ~ StLe 


1060 x 4.22 % 19.36 


; * ace 
“ 86600 Pe 
® Ley we. ¥) PER STROT 
173200 , at 
C S3F TOTAL (2 STRUTS) 
“Urt+7 
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MAIN GEAR Determine cylinder length: 


P. L= Pa (UL + 10) 
850 L = 187 (L + 10) 


L = 2.825 from static 


= .825 from compressed. 


Pressure, fully compressed 


P.L 850 x 2.825 
pel > ed — Re gat Ec el diy 
Ls 7825 
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P. = = 2915 psi 


Load at any stroke 


A. _ 2915 x .825 x 19,36 
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= PA = 2k 
LOAD = PA = =¢56 
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ma | | = 93000 - ITI-5 
S0-525 TOTAL (2 STRUTS) 
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In the simulation model all forces in the lateral plane 


ent 
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are assumed to be «ero, hence the load equations for both 
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nese and main gear are multiplied by 2 to give total load. 


(2) Damping Damping is provided by the restriction of — 
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the flow of oil through an orifice with a sharp edged entrance 
g on the lower side, and a rounded entrance (radius = .78 x plate 


thickneas) on the upper side, Damping is in both directions, 
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adding to the pneumatic force as the strut is compressed and 


subtracting as it extends. 
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The properties of the oil used in the shock strut are not- 
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known and no inforwation cencerning damping characteristics 


of the shock strut is given in the reports listed among the 
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references. For the purposes of this study, it was assumed that 


wen, 
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the oil had the following properties (at 60°F) 


4 = 500 centipoises (absolute viscosity) 
S$ = .88 (specitic gravity) 


Thus p = .85 x 62.4 = 55 #/£t2 


5 


u = 500 x 2.09 x 107° = 3072 


Slugs/ft-sec 


R = Reynolds number 


Vdop 

) a — where d, = orifice diameter 

_ | V = velocity through orifice 
| _ d, = 0.75" 

: “d, = 4.122" 


’ 
+ 


‘In terms of piston speed V’ _ 


~ {4.122\ (4.122 55 ’ 2 
. = (323 (An ) (=) — 
lo” 4 | = 

= 104 v7 


é 
V dis the compression or extension rate of the shock strut 
and as equal to the sink or rebound vate (in ft/sec) of 


é 
the aircraft. Values of V for which the damping force is 
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significant (say more than 1/10th the air compreseion force) 


are above one foot per second, hence the Reynolds number in 


all cases of interest is above 10'. In the range io! <R < 108, 
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the discharga coefficient C for a diameter ratio ao/4, # 75/4. l2e" 
0.18 is approximately constant and aqual to 0.6 (see Ref. 14). 
Solving the standard flow a juation 
: 2gaP | 
oS A ye 3 . IhI~6 


for damping force (D.F.) = AP x A. psi x in?) 
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in terms of q = VIA, (in/sec x in?) : 
© where ve is the piston or stroke velocity, in/sec,for C = 0.60 
and p = 55 #/ft? = .032 #/in® 


2 
p A, A. 2 
D.F = x in vn 
2g Cc 


2 
. 1032 x 19.36 oo 
yx 386.4 x 0.602 qs c 


a 162 (4.122\? | 2 
278 = \ .75 c 


ee 002 53022 vo" = ,067 we Vv. in in/sec) 


2 


3 c 


There is no mention of a metering pin in the available 


literature, except obliquely in Pef. 8. However, this 
extremely low value of the damping factor indicates that one 


must be present. In terms of the ratio ran of tne metering 


m 


pin diameter, 4, to the orifice diameter, aes the damping 
factor can be written as: 


S01 DR CE LINER EOI RIN EOE: 


C.F. # 10 


»2 
a Cy ve 


wheres CD is the damping constant an? r om 


mn” an 


For practical values of Ton? Cy has the following values: 


19 $2.3 | 
850 100 
975 200 ! 
.985 360 | 
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For a metering pin diameter qd. = 0.950 qe the clearance in | ae 
the orifice hole between the sides of the pin and the orifice 

edge is (1/2K.050) (.75) = (.0250) (.75) = .01875", which is 

‘glightly larger than 1/64". In the simulation, the damping 

constant (C,) was chosen to correspond to 100. Thus , for 

a sink (or rekound) speed of 10 fps, the damping force is 


D.F. = 100 x 100 = 10000# (or -10,000#). 


Both the viscosity and density af the oil are functions 


of the temperatures. The variations in viscosity should have no 


effect. Althoughthe variation is quite large, viscosity effects 

only the Reynolds number, which in turn effects the orifice 
coefficient C. But in the range of interest, C is a constant. 
Density is less affected by temperature, there bathe a drop _ 
of pev~haps 10% as the temperature changes from 0 to i00°F, 
and in view of the other uncertainties of the damping force 


calculation, this can be ignored. 


“ Usually, metering pins taper from the base to a paiatc 
near the end where the pin flares inte some sort of a bulb, 
so that maximum orifice cnnstriction occurs at aither end 
of tha atroke. This variation was ignored. 

Thus, for the present simulation, the damping force 


is assumed to be given by the expression: 


DF. = + 100 v2 111-7 


The plus sign applies when the strut in extending, since 


tha damring force is downward. 
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(3) Tire Forces In the interest of simplicity, tire 
forces are not included in this version of the simulation. 
The inclusion of the tire would add a degree of freedom to 
each landing gear. This assumption can be justified by 
noting that the work done by the tire is small in comparison 
to the work done by the shock strut. Further, the work done 
by the tire has a strong effect on the dynamics of the air~ 
craft. Since the unsprung mass begins moving at a somewhat 
later time, the real effect on landing gear forces, as far 
as the aircraft is concerned, is to delay the build-up of 
forces during the landing. However, the overall effect is 
small in cases where the unsprung mass is small and hence can 


be ignored. 


III-2 Rotor Forces, Angle of Attack aad Engine Power 


Assuming that the thrust is known, the following 
parameters can be determined: 

(a) Thrust coefficient a 
(b) Collectiva pitch o a5 and A for minimum profile drag 
(c) Torque coefficient Co and Torque Q, required 


R 
(ad) Power required Po 


(a) Power ratic (PR) - required-to-available 


{f} Rotor angle of attach a, 


(g) H = Force 
(h) Rotor Lift and dray. 
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The computations are accomplished in the following order: (2) 
a) Thrust coefficient 
Cc, = Ww (TL/¥) 


III-8 
TR (QR) 2 


b) Collective pitch setting for minimum torque 
If Eqs. 72 and 74 of Reference 13 are added, the 


resulting expression is tne net torque at the rotor 
shaft and has the form: 
= 2 2 S 
Cy K, oo + Kj6, + K,AG, + K 4A + Kg) + K, III~$ 
where the coefficients Ks are lengthy expressions which 
are defined in sub-section 3. Eq. 69 cf Ref. 12 gives 


the thrust coefficient Cy: in terms of 2 and 85? 


“4c 


: 
T + ) + 8 ) roy 
AE tangy, © Ea Me . ITI~10 4 


— = t 


a 
o. is the blade pitch at the root (collective pitch) and 


3,1 


o, is the blade twist. Since C,, is known, Eq. III-16 
may be solved for A which may be substituted in III-9, 
Th's yields a quadratic in 6° Taking the first 


derivative and setting it equal to zero gives: 


D 
oF III~il 


where Ds fate, . - - 
BEE gg t by gd Eylte g + ty oi) Sy {te y+ te yg) 


atal g $ b : : - 
“28 (sat tyg * Xs yg) - 2a (ty yt Tas * tag? 
+ pee en {é, (2 et 
aa as g,7) ~ @ (te nt te yd) 
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In these expressions 


( \ ti 5 are the Bailey coefficients, Ref. 13 
6; are the co-fficients in the section drag equation 
2 
a is the slope of the section lift curve 


and Oo. is the blade twist. 


6 antennae tae an aN (a 


Since the second derivative of Cy is positive for a given 
value of C,, Eq. IlI-llgives the value of 6, for minimum torque. 
With this value of Boe Eq. III30 may be solved for i. 


ce) Torque coefficient fo. 


The torque coefficient may be calculated by means of Eq. III~9 


oe RA aE einen tn ete ne Ar Rename Oe 


using the values of 6, obtained from Eq. IIIdl, and A from III-10. 


The coefficients of III-9 are 


Ky = (Sat, graty gt(Sat. graty git So ts a6 


Rm 
tt 


re ‘ 

2 = Sys 3 + Ayty g + (Sats graty 5) Oetzaty g 

Ky = (Oat, gu aby o) * (Sats grat, 3) TII-9a 
Bg eeoite 2: So ghesg Bhagat ] 
7 = . 2 
| Ke ™ Sots.2 * Sats. ghe * Sats 16% 
i tTerque is given by the expression 
2 2 

Q, * Cy (rR (MR}“ R) IIr-14 

Q,2 
a) Power required P, ™ yxy shp. YIr-15 


Maximum power for tho CH~53D, 4 function of temperature (°C), 
is taken from Figures 4.5 and 4.6 in Ref. 6. 


Pas » 6422 shp ot < 15sec 
= 6888 ~ 29,1t LS*C>ersseec ® 1008 nF 
= $140 bt > $9°e IYL-16 
TXI~24 
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e) The rotor angle of attack . ' 


Eq. 706 of Ref. 12 yields: 
Z h Cr 
Ge, = are tan D . + ween aKyh/2— Tri-17 


£) The H Force the profile drag-lift ratio is 


expressed as 


é § 
D = gat 6 jot | 
(Z). 2h p 3. a \te 24 bs 6,39. 75) 
gad (e, od? +t, (dO og + ty p0%D, | 
a \"6,5 6,6°°.75  “6,8°.75 


It is also equal tc 


T cos ea H sin On 


pote. 


ray 


D) H cos a. - T sin a 
ae 


Thus -T(Sin a + D ces %c) 
He c L 

cos pa + D sin @ 
«dey 
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g) Rotor lift and drac The rotor lift and drag: 


Ly ws T cos a. 7 BR ein ae III-193 


‘ Cb. vw KH Cos as T sin a. ITI~-206 


Ili-3_ The Rotor Hub Moment . 


The moment due to Elapping hinge offset ia given in 
Ref. S {(p. 49) as 


: 3 
M, = Ka)~By¢) - ($e) (=*) [) -(F) | Uxl-24 
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The term ay is the backward tilt in flapping and is 
given by the expression (Eq. 65, Ref. 12) 


= t A+ t ) 8 


a, = ty), 1,590 * *1,6% = t2,44 * t1,5°, 975 


In most forward flight conditions Q, < Bigs where Bis is the 
forward cyclic, so the first term is negative, pitching nose 
down. The term K is the hub pitch constant, ft-lbs/deg: 


eb2*Mg 
os 


The second term adds to the first if 6 is negative (pitching 


nose down), but the whole term is usually small and can be © 


ignored. Thus the hub moment is 


a : ebQ'a 

“q Ko 

4 For the CH~53D: e = 2.0 
3 b = 6 

4 5 = 184 
4 Q = 22 


Serena Re See eo 
PPB ad ANS nD IRC EI AE AEDST IE OR OIL LS 
4 


oe ee 


=a 2X §& X 22 x 184 ) 
a se 1 1s Ili~23 
= 5,35 x10° (a) ~ By) 


Thus M 


1s 


The value of a, ~ By, may be as high as ~-3° 50 My might be 


about « 238606 ft-lis. fur cha CH-<$3u. 


Tit-4  “uselage Forces and Moments 
The equations for fuselage lift, drag and pitching 
moment are based on wind tunnel data for the CH-53D prasunted 


= in Refs. 10 and 11, and shown graphically in Ref. 5. A 
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polynomial was fitted to the data for the range of fuselage tC) 
angle of attack from -1€° to +16°. The resulting equations 


are * 


LIFT Lp /q = 15.0 + 487 a, 1II-24 


42.696 - 11.45 Oy 
+ 6.423 a2 L1II-25 


DRAG sD,./q 


FOR THE RANGE -20° < a, = + 20° 
MOMENT M./q = -450 ~20 S Op = 712 III-26 
< q 
M./q = 58.8 a, “12 - a,- + 16 TII-27 
< 
Mp/q = 1000 +16 5a, +20 grreog 


The lift data is for fuselage and tail. The tail incidence 

equals 3°. The drag data is for the complete aircraft, in- Oj 
cluding rotor head, antenna, cocling losses, etc. The moment 

data is for the complete aircraft, minus the tail. The full- 

scale data shown in Ref. 5 has been transferred to a c.g, 


at station 386 and a waterline station of 161.4. 


IIx-S Tail Forces and Moments 


The horizontal tall lift esefficiant is 


“i =a (A, - a.) 


TiI-29.1 


where the angle of attack tp ds 
an ™ "+ + ib ~ , + 7 


and % is the garo lift angle. The slope of the lift curve 


ITI-29.2 


ae#cC 
bat “3 


1+ VE) : +i 
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where a, is the lift curve slope for the airfoil section 


and is taken to equal 5.73. For ah aspect ratio equal to Z 


2.5, a = 2.88/RAD or 0.05/DEG. : 


Op. the angle of attack of the tail, is equal to 
the fuselage pitch angle Ser the tail plane incidence angle 
ine the downwash interference angle Eq and an additional 


angle of attack change produced by the fuselage pitching 
motion and equal to 


8 ¢Xy 


where fT is the distance of the tail to the c.g. The angle 


€, is discussed in the next section. 
| e 
S| SUN ee, OO 8 iTI-29.4 
EO Mop Lp 
4 4 Tail drag is a part of fuselage drag and included in 
4 Equation 111-25, 

IL-6 Taterference Angles 
4 The generalized equakions ef motion santain the 

= .dewnwreh intexrfarence angles | 

S " | - (i) fens ~Soter on fuselage 131-30 


(2} tn Reker on tail TIF+3) 
which apply in this cage, aa wvel\ as several others (Cn, 
Ege Eye Epps and Song! which do nok apply. The formulas 
feyg ONG ty SYO taken from Ref. §, section §, 

fpye % Tan a, 
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Some Consideration was given to the use of the formula 


“ -A 
where Kop igs a function of the rotor wake angle 


x= a, + tan”? E u | 


The functional relationship is shown in Ref. 4, Section 5, 


Figure 1, For a range of x from 40° to 80°, Kor {and Ren) 
ranges from 0.37 to 1.46. Since no information was availa- 
ble on the wake angle x for the flight conditions of interest, 


it was decided to defer this question, 


IfI-7 Automatic Flight Control System (AFCS) 


The simulation model of the AFCS is mathematically 


identical, except in one minor respect, to the AFCS instalird { , 


in the CH-53D aircraft. The exception is due to the finite 
difference nature of the AFCS simulation which @elays the 
feedback one time step. The pitch transfer function ef the 
APCS feedback locp is derived as follows: 
By, Ryle + Kite 

where By, is the cyclic pitch of the swash plate 

* ig the pitch angle ef the fusolage and K, and Ky 
are the rate and proportionel ssine. The transfer function is 

ait “x, (rs eo 


x F 
R 
where T Ae). The feedback gains for optimum performance 


hate been Jerived in Reference $, for th- aircrafe Sys tam 


go tae + 2} "cyclic/* fuselage 
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-and the value of Bis is calculated from the values of 8 
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The question of time delay referred to above is of 


interest. If the APCS simulaticn were described by con~ . 
tinuous (analytic) differential equations, the time delay 


ae and its 


Gould be represented by the transfer function e- 
effect could be readily determined by the various techniques 
cf control system analysis, say by a Nyquist plot. However, 
when the AFCS simulation employs difference equations, the 
control feedback is delayed, i.e., 

(By ,) = KX, (Bp) + K_ (85) 


P 
tat, +T t. tet, 


£ 
and 6, at the previous time step. 


Figure ITI-2 shows the effect of step size. The curves 
trace the motion in the x-z plane. At the end of seven 


intervals (8.4 seconds) the coordinates are: 


fh Sih 2. wh 
010 849 +20 3.85 
£625 336 .70 «3,00 
63) 168 1.02 3,38 
.100* ad 2.30 4,10 


150° 69 4.96 6.90 


tcurves not shown in Figura Iif-2, | 
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It can be seen from the results above and from Figure 
III-2 that the error due to the finite difference representa- 
tion of the AFCS is not large if the step size is less than 
050, for the hover conditions only. Conservatively assuming 
that all the differences in Figure III-2 were due to the 
delayed feedback error, then the error incurred in increasing 


the step size from .010 to .025 would be « 


c= 2(+70-.20)* + (3.85-3.00)7 _ .987 
. 336 336 
= 00294 


which is less than h but greater than nh’, It can be shown 
by the examination of a simplified version of the equations 


of motion that the delayed feedback error is of order Nh, but 


“not wh’, where N is a positive number which is a function 


of the x force. 


TiI-8 The Ground Effect 


cele 
Figure III-3 shows the variation of the ground effect 

parameter 

A = rotor thrust in ground effect 

rotor thrust out of ground effect 

as a function of the ratio 2/R 

é = distance of gs.0r head above ground 

Ro rotor radius 
in torms of the ratic V/U, 


Ve forward speed (~t1)) 


fom 4/2 
Uo ( a « induced velocity in hover 
20 oR? Be 


Va tip loss factor # 0.9? 
(Tie2es 
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Figure IlI-3 Ground Effect Parameter A vs. Z/R 
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For T = W = 33500 


33500 os 
U5 a 2 2 
27) (0.00238) (36) “(6.97) 


= 76.7 fps 


Since the launch/recovery operations are limited to an initial 
or terminal speed of about 25 knots ship velocity and 45 
knots wind velocity (rhe maximum allowed oe CH-53D wind-up) 
V = 70 knots at take-off 
= 118 fps. 
and V/8, = 1,55, 

The parameter A gives the reduction in thrust due to 
ground effect. Ground effect will alter each of the variables 
computed above, since the effective thrust is reduced. Since 
the LARC~I calculation begins with a value of thrust (in 
terms of T./W), and this value ir multiplied by A {as a 


function of F/U,) then the values for Coe C Ons ong 


Q,' “a Cc 

R 

H, Lp and D.. all include the ground effect. However, it should 
be notec that the unbalanced moment due to a rotor being only 
partially in the ground effect, as when the rotor extends over 


the edge of the dack, is net included. 
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IV. GENERAL DESCRIPTION OF THE LARC-I SIMULATION iS) 
The LARC-I (LAunch and Recovery Capability) simulation - 

provides the necessary coupling between the motions of the 

ship, discussed in Section II, and the aircraft motions, 

discussed in Section III. The motions of the ship are con- 

fined to the x-z plane (i.e., to heave, pitch and longi- 

tudinal velocity) and those of the aircraft to longitudinal 

and vertical velocity and to pitching about the aircraft 

center of gravity. 


Physical interaction of the ship and the aircraft occurs 


Only through the landing gear forces and moment: 


2 =z- (F Fee) 
LG LG, LG,, 
Xig = brake force = + up 24 a 
= F X TPS 
PUG LG, “LG, LG “LH, 
(up = coefficient of friction, usually 0.3) 
Pug and PLGy are, respectively, functions of Sy and Sq! the 


extension or stroke of the nose and main landing gear shock 
strut cylinders which are in turn functions of the relative 
aireraft and ship positions. The equations for stroke 
are shown in Appendix C, page C-1. It is only in these 
equations, in the torms ie (vertical position of the ship 


two aircraft lLauneh/landing points) and . (velocity of 
LP 


the same points) that the ship motion interacts with the 


aircraft and the interaction through the landing gear forces 
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persists only as long as the aircraft is in contact with 


the deck. 


The variables 2. and Z are in turn harmonic 
“Lp itp 
functions of time which include the coupled heaving and 


pitching of the ship in response to waves of a random sea. 


The sea state is, in this simulation, characterized by 


a single parameter, the significant wave height his3 which 


may take values (crest to trough) of up to 100 feet in sea 
state 9 conditions. 
Further coupling between the ship and the aircraft 
occurs indirectly through their relative motions. Since 
both the ship and aircraft motions are referred to a body axis 
system it is necessary tc refer tne motion of each vehicle 
to known fixed inertial axes (see Figure IV-1) in order 
to obtain the re‘ative motion between the vehicles. The 


n2ecersary transformation equations are given in Appendix a. 


Iv-\_ Structure of the LARC-I Simulation Program 
Figure IV-2 is a simplified block diagram of the 


LARC~I simulation program, The simulation may be described 
as 4 fixed control non-steady state model, Each computer 
run simulates a scqment of length T (seconds) during which 


the program computes the motion at N «= 7 
b¥ 


instants (aT is the time step) ducing which 211 input 


vatiablos including the controls remain fixed. If all air- 
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are set to zero, the aircraft is placed on the deck at a 


distance J from the ship's center of gravity. 
LP 


The sequence of operations for the simulation which 
uses a fourth-order Runge Kutta scheme to numerically inte-. 
grate the equations of motion is shown in Figure IV~2. At 
each time instant qT. the force and moment modules (2), (3) 
and (7) compute the forces and moments based on the attitude 
and velocity of the aircraft at the end of the previous 
time step, T.-At, and, if the aircraft is in contact with 


the deck, the position and velocity of the ship, block (7). 


At each time Ty. a value of Be is calculated at block 9 
Ss 


based on the value of a and ie at the previous time TAT. 


If the thrust/weight parameter (T/W) is greater than 
one, the aircraft will ascend, if less than one, the aircraft 
wil], descend. The thrust may be delayed Dy seconds by assign- 


ing a value to the delay parameter D This allows the take~ 


D* 
off sequence to start at some specified point in the ship 


motion cycle. Experience has shown that the parameter Dy, 


is less important in the landing sequence. Also, the thrust 


may be specified to increasa linearly from zere to (T/W)W 


pounds in Dy seconds for use primarily in the take-off 
sequence, 
The initial longitudinal cyclic, Be may be specified, 


nd 
$ 
However, tha Automatic Plight Control Sy8tem will continuously 


Vary By according to the equation 
“S$ 
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Q Bi, = K.8, + Ki, (8—-8,) 


where K. is the proportional gain, K, the rate gain, O- 


R 
the fuselage pitch angle, and B ; is the swash plate tilt 
“Ss 


angle relative to the rotor shaft. The effect of assigning 
a value to Bf may reduce the magnitude of a transcient 
S 


motion, since otherwise By is assumed zero. 
s 


fo) 

The trim pitch angle is the fuselage pitch corresponding 
to near equilibrium conditions at a thrust equal to (T/W)W 
and for all other vaiues of the initial conditions. For a 
given gross weight and c.g. position, 0, is a function of 
forward velocity. If a is not known for a given velocity, 
it may be found by a trial and error procedure. Since the 
LARC program outputs all forces as functions cf time, the 
trim angle for a given speed may he easily approximated using 
the data from a few runs. Although trim pitch is specified, 
the aircraft in the simulation is rarely, if ever, in a ~*ate 
of trim, or steady state flight, although the changes about 


the trim value will net be large. The thrust program is shown 


- : 
| 
4 
BS 
Fg 
; 
k 
: 


in binck 4. The integration scheme is shown in biock 8. The 
dashed Glock 11] refers to tha model for gusts or turbulence, 
which, although very siinple, was not included in the present 


version of the simulation. 


W-2  Qperation of the LARC-I Program 


Tha following remarks are intended to jilustrate the 


operation of the simulation. They are not intended to 
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serve aS a user manual. Operation of the program is quite 
simple. A user with no prior knowledge of the program itself, 
or prior background in programming, can make full independent 
use of the program after only very brief instructions. 
The program is designed for interactive operation on 
a time-shared computer using a remote terminal. The user 
communicates by means of the terminal keyboard. The program 
response is by terminal print-out or graphic display. Equip- 
ment used in this project was a Honeywell H-1648 time-shared 
computer, accessed by various remote terminals such as 
Tektronix 4010 and Typagraph 3 and DP-30. 
The LARC-I program written almost entirely in FORTRAN 

IV, is essentially machine-independent. Certain output, 
plotting and file processing routines which are required by 
the particular interactive graphics terminal (Tektronix 4010) 
are features of the H-1648 time-shared operating system. 
Since these features are common to most time-shared systems, 
operation of the program or computer systems ether than the 
H-1648 can be accomplished with relatively small changes in 
the coding, Machine independency wae a major goal of this 

reject. 

| All processing on the time-shared system is done using 
files whish contain programs or data; each file is created 
independently. The files are accussud by reference to file 
namo, either internally in the program or by means of an 


interactive command. Data files may be created in thair 
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entirety, may be searched for specific elements, or may be 

changed hy altering specified elements or groups of elements. 
Subsequent to the creation of the data files referred 

to below, and to the compilation of the LARC-1 program, 

the sequence of interactive commands may appear as follows 


(the user replies are underlined): 


TYPE AIRCRAFT TYPE *** CH53 

The user response leads a file entitled “CH53", containing 
the input data shown in Appendix C. The user may select 
any data file which he has previously created, each repre- 
senting a particular a‘treraft configuration. tlements af 
this file may not be altered during run time. 

tt D 

The user response loads the parameter file DEFALT. The 
values of this file may be changed during run-time in the 


Manner described below 


wt 


qe- 


The user resporrs "I" for INPUT, signals the intention af 
changing Values in the parameter file DEFALT. If this re- 
spouse ia not made, the OEFALT file is loaded ag compilad 
prior ta fun<tine, 

PARAMETER CODES « YES UR NQ | NG 

tf the user response ia "YES", ihe parameter aude are printed 
as shows in Table IV-r, These oodes ageign a nuchor to be 
used ia the roponses tiated below, 


INPUT COME AND VALVE 21,0.2% 
we a Ata 


: Bio Sac chee weeks 
AOR engi s en dats pel Satna Sse tha: SS 
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Af new PARA 4 a 


The user has assigned a value 0.05 to parameter 21 (the 
run-time). The user may assign values to any or all of the 


parameters listed in Tabie IV-1 in any desired order. 


INPUT CODE AND VALUE 25 
This response is required to terminate the input sequence 
following the command "I" and return the program to the 


command mode. 


ke OT 


The user response is "T° for take-off but may be "L" for 
landing. This initiates the computational sequence which 
continues uninter:‘upted for N = T/AT steps. Both T (run- 


time) and AT (time step interval) are parameters and may 


be altered for each run. 


ner 


ex P ef 
‘In the present version, the user may command outout in 

tabular form, "P“ or graphical form, "G". 

ROW MANY CURVES (1 - 3) 3 


1. ~ CURVE @ ** 


im 


2. - CURVE 


at 

* 

= 
~~ In 
in 


3, = CURVE & ** 


{ 


Hore, the usor specifies the curves to be printed or dis- 
Played graphically in any de.ired order, in groupa of one 

to throes curves. Apnerdix a santains exatples of the printed 
eutput. Table IV-2 lists the output paxametera whieh may 

be priated on seaviand. These 23 files are jloaded in their 


entiroty each time the program executes through the cosputation 
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sequence. The user may limit the output in any desired manner 
by specifying from 1 to 30 curves in groups of one to three. 
After the output of each group, the program will print the 
command signal "**" at which time the user may respond 
with: 
Por G continuing the svutput sequences as described above 
DB initiating another run in which parameter 
changes may be made as described above. 

Other commands which might be made are: 


N for NEW CASE which re-initiates the input 


Sera: eer meer enter etree ananaa eee S 


cycle from th-: beginning. The program response 


to this is: 
TYPE AIRCRAFT TYPE 
allowing the user to load the file pertaining : 
to another configurat) ix model of th. ame 


aireraft, ar a different aircraft. 


~“ 


L for LANDING which sets certain program switches 


pertinent to the lending calculation. 


v3 


far ENG which stops execntion. 


Regardless of the number of steps Nos T/AT, the program 

stores only 50 values. Al} 506 values ara used in graphing 

the variables ghuwn .n Table IV-2, but only 25 valves ara 

wrinted. if © >aer dezi: * more accurate rerolution, he 7 
must #2: Tish this by ccareasing the nutoer of steps 

foe T/AT, by shortening the run time T. N must equal at 


least $9. i 
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Individual runs may be continuations of a prior run. 
Teble IV-3 and Figure IV-3 illustrate; the procedure. The 
values in Table IV-3 for cas2 I are taken from Appendix A, 
Series II, case 1. The run time for this case was 15 
seconds. The analyst chose to change thrust and fuselage 
trim at t = 9.6 seconds. The initial conditons for case 1.1 
are those which prevailed at t = 9.6 in case 1, plus the 
control changes. ‘his procedure may be repeated any 
number of times to compute the flight path of a maneuver, 


during each segment of which the controls are fixed. 
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TABLE JV-1 
PARAMETER CODES 


AIRCRAFT INITIAL CONDITIONS 


1. A/C INERTIAL X POSITION (x) 

2. A/C INERTIAL Z POSITION (z) | 
3. A/C PITCH RATE (8) : 
4, A/C PITCH ATTITUDE (8) ‘ 
5. A/C VELOCITY - X DIRECTION-BODY AXES (u) ; 
6. A/C VELOCITY - 2 DIRECTION-BODY AXES (Ww) | 


AIRCRAFT CONTROLS 


7. THRUST TO WEIGET RATXO (T/W) 
8. CYCLIC PITCH-LONGITUDINAL (B, ,) 
9. FUSELAGE PITCH TRIM ANGLE (9,) 


SHIP MOTIONS 


10. SHIP HEAVE AMPLITUDE (ny) 
ll, SHIP PITCH AMPLITUDE (ng) 
12. SHIP ENCOUNTER FREQUENCY (u) 
13. PHASE ANGLE - PITCH TO HEAVE (a) 


14, INITIAL SHIP X CG INERTIAL 
15. INITIAL SHIP %4 CG INERTIAL 
16. INITIAL LAUNCH/RECOVERY FT ~ X 


1?, INITIAL LAUNCH/RECOVERY FT ~ & 
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TABLE IV-1 (Cont'd) 


SHIP VELOCITY (v_) 


WIND AND SEA STATE CONDITIONS 


18. WIND VELOCITY (Vi) 
19. SIGNIFICANT WAVE HEIGHT (hy 73) 
20. INTEGRATION CONSTANTS 

21. RUN TIME (T) 


22. TIME STEP (AP) 


PARAMETER CODES 
THRUST TIME CONSTANTS 
23. THRUST DELAY TIME (D 
24. THRUST BUILD-UP TIME (D 
PROGRAM CODES 
25. TIME AT INITIAL STEP (T )} 


26. CONTINUE 


iv-}2 


OE ee RE ee ee Ee a ae ee ee a ae a ae ee eae eg ERE NR ee ne ee ee aR hy eae, NE Ee eR ORO OS aT OR, a ge OR ROL ee EE ye 


TABLE IV=-2 
OUTPUT PARAMETER LIST 
CURVE CURVE 
NUMBER LABEL PARAMETER 
1 U A/c “A/C FLIGHT VELOCITY BODY AXES 
2 W A/C A/C VERTICAL VELOCITY BODY AXES 
3 THETA DOT A/C PITCH RATE BODY AXES 
4 THETA A/C PITCH ATTITUDE BODY AXES 
5 X(I) A/C A/c X POSITION INERTIAL AXES 
: 6 z(I) A/C A/C Z POSITION INERTIAL AXES 
7 XISCG SHIP CG X POSITION INERTIAL AXES 
g ZISCG SHIP CG Z POSITION INERTIAL AXES 
9 FWD CP NOSE WHEEL DECK CONTACT POINT INERTIAL AXES 
10 AFT CP MAIN WHEEL DECK CONTACT POINT INERTIAL AXES 
a2 FLGN A/C NOSE GEAR LOAD LBS 
: 12 FLGM A/C MAIN GEAR LOAD LBS 
13 BIF CYCLIC PITCH-LONGITUDINAL RAD 
14 ALPHA POTOR ANGLE OF ATTACK RAD 
15 GAMMA CLIMB (9R DESCENT) ANGLE RAD 
16 THRUST THRUST LBS 
17 aE ROTOR 2 FORCE (BODY AXES) LBS 
18 XP ROTOK X FORCE (BODY AXES) LBS 
19 ZLG LANDING GEAR FORCE LBS 
20 XLG LANDING GEAR SRAKE FORCES LBS 
22 MLG LANDING PITCHING MOMENT FT-LBs 
IVv=14 
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CURVE 
NUMBER 
22 

23 

24 

25 

26 

27 

28 


29 


CURVE 
LABEL 


MGL 
NGL 
LF 
DF 
LFUS 
DFUS 
LT 


DT 


TABLE IV-2 (Cont'd) 


PARAMETER 


MAIN GEAR STROKE 
NOSE GEAR STROKE 
LIFT ROTOR 
DRAG-ROTOR 
LIFT-FUSELAGE 
DRAG-FUSELAGE 
LIFT-TAI i. 


DRAG-TAIL 


() 


INCHES 
INCHES 
LBS 
LBS 
LBS 
LES 


LBS 
LBS 
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V.._ RESULTS 

Over 350 computer runs were made, 34 of which are shown 
in the tabulations of Appendix A. The principal object of 
these runs was to test the program and to gain éxperience 
in its use. Thus, these runs do not necessarily represent 
a systematic evaluation of the aircraft/ship compatability 
problem. Budget constraints prevented an in-depth parametric 
survey of all aspects of this problem. In fact, the “ship” 
itself is only an approximation, since no detailed and 
accurate information was available to model a specific ship 
design. However, important conclusions have been drawn 
about the effects of ship motions on aircraft control and 
stability: 

(1) Even for the most severe ship motions (a wave 
height of 30 f£°., an encounter frequency of 0.3 rad/sec., 
Series III, Case 11) the CH-53D control margin is adequate 
for norma] take-off. For ship cncounter frequencies from 
0.3 to 0.98 rad/sec., corresponding roughly to periods 
of 2.695 to 9.0 saconds, the ship motions arc quite slow 
in comparison ta the aireruft motions. Series III, case 4, 
Yopresents a severe ship metion. In this cese, the ship 
is heaving in a 30<ft. wave with an amplitude of 18.3 feet, 
and a pericd of 3 seconds. The ship is at garo amplitude 
att « o, The maximum aircraft rotor thrust is attained 
at te $ seconds and the alrecraf: is aff the deck by t © $.0 


goconds. Desrite the large ship oexcurs.ons, there is no 
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excessive build-up in landing gear forces. The aircraft 
pitches forward about 3.5 degrees during the first three 
seconds, lifting the main gear off the deck. The nose gear 
force remains low, however, and vanishes at about t = 5.0 
seconds. At no time during this 5-second period does the 
pitching acceleration of the aircraft exceed + 0.105 rad/sec. 
(2) In ail the runs tabulated in Appendix III. take-off 
occurred from a rising deck, at a time when the deck upward 
velocity was near a maximum. Other runs, monitored visually, 
were made for take-off from descending decks, also at 
maximum velocity. In no case were large or destabilizing 
moments transmitted to the aircraft, Case 1 of Series III 
is typical. Here, the deck is ascending at a rate of 
# (15) (.68) 2= 7.0 ft/sec at t = 0. It is assumed 


2 
1/3" i) 
that the aircraft is released at t = 0 at a distance of 7.5 


/2 h 


inches apove the static gear «xtension. Even so, the maximum 
total gear load is only 46500 por-ds at t = 0.6 sec, {A/C 
gross weight is 23500 lbs.), and the landing gear pitching 
moment ig -104000 in.-lbs. nose down. The pitching velocity 
(THETA BOT) remains very low. 
(3) Sanee the ship piteh amplitude was small (8.492% or 

. less) all runs were made with take-off from a point above . 
the ship center of gravity. In the worst case, the ship's 
pitching motion, leading the heave by as little as = radians, 
would add luss than 1.3 feot to the Cotal vertical aneuaes 


ar either bow or stern, 
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(4) The gains for the AFCS pitch channel were determined 
for the CH-53D and agree with the optimum values calculated 


by Sikorsky (Ref. 5). These values are: 


K, = .4 rate given 


Kp = .8 proportional gain 


(5) In its present form, the pregram can be used to 
determine gains, stability margins or trim conditions, but 
the procedures are not “optimun", Experience has pointed 
out changes which can be made to effect rapid conyergenge 
to optimum values by use of acaptive methods. Once optimum 
gains and trim conditions are identified, stability margins 
can be readily identified by an analysis of the motion. 

(6) Although roll and yew are suppressed in this model, 
the simulation yiclds considerable data which provides in- 
sight into the roll hehavior. The significant variations 
in 4 (the inflow ratio) and By, (the longitudisal cyclic) 
at frequencies in the neighborhood of 2 to 3 rad/sec., be- 
cause of the effect on side foree, conld induce critical 
eapupling for shin motions in the sama frequency range. The 


coupling effects stem fren the presence of \ in the side 


foree equation and in the third harmonic flapping term b,. 
a. 
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APPENDIX B 
GENERALIZED AIKCRAFT EQUATTONS OF 
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GENERALIZED EQUATION 
FOR 
SIX DEGREES OF FREEDOM 
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APPENDIX C 
SUMMARY OF 
CHARACTERISTICS, FORCES AND MOMENTS 


OF THE 
CH-53 HELICOPTER 
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SUMMARY OF 


CH~53 CHARACTERISTICS AND PARAMETERS 


CHARACTERISTICS 

Rotor 

R RADIUS, FT 36 
b NUMBER OF BI.ADES - = 
g ROTATIONAL SPEED, RAD/SEC (NOTE 1) —21.7 
6 SOLIDITY RATIO | 1150 
e FLAPPING HINGE OFFSET, FT. 2.0 
e, BLADE TWIST, DEG. ~6° 
be ROTOR SHAFT INCZDENCE, DEG. (LONG.) 

(NOTE 2) 5° 
Y BLADE MASS FACTOR (LOCKE NO.) 12.55 
M. FIRST MASS MOMENT OF BLADES, 

SLUG ~ FT 183.6 

lie x,y COORDINATES OF 7.50 
1 ROTOR HEAD, FY 1533 
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CHARACTERISTICS 


TAIL 
2 
A, AREA, FT 40 
i INCIDENCE, DEG. 2.0 
& ASPECT RATIO 2.5 
1, ” ~9.85 
COORDINATES, TAIL AERODYNAMIC 
a CENTER 56086 
LANDING GEAR 
1 
z 8.16 
i g n : 
i NOSE 
“Fy COORDINATES OF GROUND : 
CONTACT POINT 
1 1 
Zy 8.16 
1 “Sy MAIN 
x 1 20.0 
Fm 
AIRCRAFT PARAMETERS. 
W WEIGHT, LBS 33500 
rt MOMENT OF INERTIA IN PITCH 
YY SLUG ~ Fr? neous 
e AM@IENT TEMPERATURE, °F 59.2 
Ko. @ ENGINE TORQUE 126. 
AFCS CONSTANTS 
K. PROPORTIONAL GAIN 0.8 
Kp RATE GAIN 0.4 
Ke SUMMATION GAIN 0.0 
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(4). 


NOTES 


Rotational speed, power on, is maintained between 21.8 
rad/sec., max, and 19.4 rad/sec. min. (Ref. 5). 
Inclined forward. 

Positive cyclic tilts the swash plate counterclockwise, 
Parameters can be changed interactively at run time. 
The procedure is discussed in Chapter III. The 

values in the table are used unless specifically 


changed at the beginning of each run. 
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SUMMARY OF FORCES AND MOMENTS 


LANDING GEAR (Section III-1) 


NOSE GEAR LOAD VS STROKE (S.) 


LG, 


173200/(S,, + 4.22) IiI- 4 


MAIN GEAR LOAD VS STROKE (5S 


Ly 
LG. 


DVS STROKE (5,4) 


93000/(S,, + 0.825) III-5 
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PT ITY AS _A FUNCTION OF A/C_AND SHIP VELOCITY 
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stroke velocity. (Subscripts: N=nose, M=main) 


FUNCTION OF AIRCRAFT/SHIP POSITION 


gear ¢'s are perpendicular to A/C 


x-axis) 
- Lie ~ Za + 1]- 8 
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ac LP. LG, 
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DO 0 ae ee 
A/C 


Distance from cg to landing gear extended, 
parallel to z axis, same for both gears. 


Launch/landing pad contuct point (inertial frame) 


Distance in x direction, landing gear 
eontact point te eq. 


Pitch angle, fuselage, positive nose up. 
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ROTOR FORCES, ANGLE OF ATTACK, ENGINE POWER (Section IiI-2) fe 
1, THRUST COEFFICIENT 


TR p (QR) 2 


CT = Thrust coefficient 


W 


i) 


Gross wt., lbs 
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I, = Rotor radius, ft. 
R 


Rotor radius, ft. 


it 


¢ = Mass density of air, slugs/£t° 


A : Q = Rotor rotational velocity, rad/sec* 
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4. POWER REQUIRED 
~— eo a 
Pp = (Qp%) /550 III-1: 


Power available (maximum) 


Pa = 6422 shp. t < 15°%c 
max 
= 6858 -(29.1)t 15° >t > 59°%c III~16 
= 5140 t > 59° 


s. ROTOR ANGLE OF ATTACK 
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For these studies, the following values were used 


é # O.0087 


. ° 

1 4) a =G. 9216 

3 $4 « +9, 400 

: & is the siopo «£ the section lift curve = §,73 
q o. is the blade twist 
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Rotor Forces (Continued) 


3. TORQUE COEFFICIENT fo ae 
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4 where o = solidity ratio (3 
4 } = inflow ratio [({V sin a-v)/QR 
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ROTOR HUB MOMENT Section III-3 
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ay = to + ty 5% + ty 62 
where e = flapping hinge offset, ft. 
Bis = eyclic pitch, rad. 
a= section lift slope 
M., = Mass inoiment of blades, slug - ft. 
a, * second Fourler flapping coefficient, (tip- 
path tilt} in flapping equation 
B= a, ~ ajcos ¢ ~ bysin @ > s> 
(Ref. 32) 
FUSELAGE FORCES AND MOMENT (Section Iiz~4) 
DRAG Op/q = 41,696 = 21.dS ay + 4.423 of IY1-25 
for -20° © Sp 5 #20° 
; tp F Uy fpug “ YO® Piselaga angle of attack 
feue * interference of retor on fuaulaga 
* 
qe dynanic pressure « b/gsy 
¥ = angls of climb 
# are €an ts 
i=9 : 
islet aie mwa ary. noe Reeth £. 0 es wae kok ce TR Ee IRE NE ee LP TERE CO Oa NE een OPE PE RTT TIRE | eRe UIP EL Nee MNT VT Te Corner Y St, sy 


ap sis ap taed Seb cpap cia et 


s Teas Sih ee ta err tre a rena. 
beara gS fe PAN MRA OR. emit SU La cemacbeden inh Lindl ies amber MS LAL DSi ee 


Fuselage Forces and Moments (Continued) 


u = velocity in x iirection > 


w= velocity in z direction : 
v= f(a? 4 wy . 


Moment M,/q = ~450 ~20°< Gp< ~12° TII-~26 


=58.8a,, “l2 < os < +16° ~TTT~97 
TII-28 
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TAIL FORCES AND MOMENTS (fection III-5) 
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AFPENDIX D 


THE REPRESENTATION OF RANDOM SEA 
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D-1 THE EQUATIONS OF SHIP MOTION 


If it is assumed that the motions of an arbitrary ship 
are linear and harmonic, the six linear coupled differential 
equations of motion can be written in the following abbreviated 


form: 


where the six displacements are each denoted by n, where the 


subscripts from 1 to 6 refer respectively to surge, sway, 


heave, roll, pitch and yaw. 
It is not our intention to discuss the derivation or 


solution of equation D-1 or the meaning of the coefficients 


M, B and C. Experience with these types of equations are. 
assumed and it is sufficient to say that M is a matrix cf 


mass terms augmented in sume fashion by the added mass of 


the displaced water, B is a matrix of damping constants and 

Cc is a matrix ef hydrostatic restoring sosfficlents, Far 

& ship with lateral symmetry under the assumption noted above, 
the surge, heave and sway equations may be solved independently 


oF the other three, and the solutions of ali six equations will 


be uimple harmenic functions 


ny a A, cos {ut + ey) ged. é De2 
“2 , i2 the forcing function, or the equatior of the wave (F ef Wt, 
" nan be expressed in the simple narmoniec furs 
a F cos wt * b-3 
pep -*% 
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where F is the amplitude of the wave 
w is the encounter frequency 


A, is amplitude of the ship response. 


We will only be concerned with equations D-2 and D-3 and for 

all other matters, the reader is referred to References 1] - 3 

and especially Reference 15 which is a thorough and exhaustive 

discussion of the derivation and solution o£ Equation D-l. 
Equation D-2 is the sinusoidal response of the ship at 

a given frequency » to a sinusoidal wave of form D-3. In the 


next section, itt will be shown that the real sea profile can 


_in no realistic way be represented or approximated by the 


form D-3. The teal sea can in no way known be represented 
either as instantaeous profiles 

2 = F(x) 

Z = F(x,y) 
or as time histories 

2 = F(x,t) 

2 = F(x, y,t) 
where the representation F is a sum, finite or infinite, of 
periodic functions, sinusoid or other, or polyncmials suci 
that the error 

em F - G(x,y,t) 
betwean the representational function F and a measured profilla 
of Wave amplituces G is cither less than a preassigned value, 
or ig a mathematical function, analytic or not, of any of 
the variables. This ds, perhaps, the brcadest statement that 


can bo made in mathematical terms of what apsears to be a truth 


D-2 
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which can be stated in less mathematical terms as follows: 


ae 
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it does not appear to be possible with the mathematical tools van! 
we have at hand to make a mathematical deterministic model 
of the profile of the sea surface which is a function of 
no more than a few variables, say 5 or 6. If the variables 
are limited to x, y, z and t, it is cieariy impossible, 
Even if othe: variables are added such as wind strength, wind 
duration, wind direction, fetch length, boundary values of | 
coasts or water depths, etc., a deterministic model is 
still impossible. 

If this is true, the engineer has a problem. In nearly 
all important cases, in dynamics, the equations which he 
deals with are equations of motion or are derived from 
equations of motion, and have the form of equation D-i. 
The statements made in the prior paragraph mean that the 
expression Ever cannot be used to represent real waves. 
The statements seem to go much further: since there is no 
deterministic way (that is a deterministic function of a few 
variables) of realistically representing a forcing tunction 
in equation Del, then there ig no way of comparing a solution 
of equation D-l] with what may occur in reality. For example, 


we can calculate the motions of a ship in the sinusoidal 


wave a. cos wtirwe can measure the motion in a towing tank 


where such a "rure" wave shape can be generated, and compare 


tha results. This is done in Reference } through 3. But 


We cannot coirpare the motion in a real sea with the calculated 


3 motion in “pure“ waves. It was said above that the statements 
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in the prior paragraph "seem" to be saying this, but this, 
fortunately, is not altogether true. There is some truth 
but it is not the complete truth. In‘what follows, it will 


be shown that a statistical (as opposed ty» deterministic) 
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use of the solution of Equation D-1 can be derived using a 
statistical representation of the real sea. The theory 

of the statistical representation of the random sea is rather 
simple. There are some serious problems in using the theory, 
beyond the scope of this brief survey, but which are largely 


concerned with how we measure wave amplitudes. It is important 


to remember that the stutistical model is completely empirical. 
The deterministic model is partly empirical. | 

We have dwelt upon this, at some length because the re- 
conciliation of the deterministic model, exemplified by 
Equation D-1, of ship motions, and the statistical model of 
wave motions discussed below, although not difficult to under- 
Stand, can lead to serious misconceptions as to what the results 
really mean. In most cases, the true meanings are probabilistic: 
“an amplitude of a given value will be exceeded only ina 
given time.“ What, then, is the meaning of choosing one 
ampla aude rather than another? Ansvers to this and other 


questions will be attempted in the foliewing discussion. 


Leh ane aot ata ie bt a RE 


A re IY, a ee a Be 


eS, ange he 
ROT Cae pasinicsnl 


os nt a ee ~ seeks at - De: a = rv Dania cekat x =. 2 Hatter ora aa ne tee 
witdhy SNS cekatee bE et ot 
+ etek A eee en a Rad TT ASP RA i eT a: 
atnnacr ene mena ree ete aA a LEAS nrrmee anes = 


D-2 the Nature of the Irregular Sea 


A typical contour plot of the sea is shown in Figure D-l. 


a _ The smallest circles locate the highest (or lowest) points. 
x One should be impressed by the disturbing irregularity of 


the pattern, disturbing, that is, to an applied mathematician. 


Figure D-l Typicul Contour Plot of the Sca 


The pieture is very Messy ao mutter which way you hang it. 

It the pattern appears to have a prevailing motien (say, 

from left te right) we might expect the crests (and troughs) 
to extend indefinitely aft right angles to the motion, They 
obvicusly da not. This attribure of the wave patturn is 
Known as short crestedness, since no one crest extends very 
far and gnethor appears not far away. Those short crests are 
nob aligned with sach other, or of the same length ov height. 


{Campare Creets a-a, b-b, and c-c}. The contours of Figure 
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-D-1 are contemporaneous and if one were to record a multitude 


of such plots at succeeding instants of time and if it were 


possible to follow the history of each contour, it would 
appear as if the various crests were meving at different 
velocities and in different directions, although the entire 
pattern 'seems' to have a wedualttag speed and direction. 
Also, crests are not very durable, They appear and disappear 
and their total life may not be very much longer than the 
apparent “period" of the wave. 
If a cut is taken along some line AA, the resulting 


profile might look like Figure D-2. 
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Figure D-2 Typical Wave Profile 


This is the same kind of profile which might be measured 
by a wave height measuring device over a period of time. 
That is, the spatial profile of a wave pattern looks similar 
to that measured at a peint (say 4) as the waves move past. 
The profile ef Figure D-2 seems easy enough to undaer- 
stand. Lot us say it was made by observing the height of 
the water on a staff fixed in some fashion to the ocean floor. 
Neediess to say, this is not how the record was obtained but 
the remarks made below apply in any: case, 
Although Figura D-2 looks simple, it is enormously 


difficult to interpret. Tha “waves" are coming from many 
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directions although there seems to be some kind of prevailing 
direction. The crest lengths are all different. Point 1 
moving past B is only the low end of a crest, but Point 2 

is a high point. Furthermore, the crests are moving at 
different speeds although they give the appearance of moving 
in groups. A stone, for example, thrown in a still pond 
induces a wave train that travels toward the shore at a steady 
"group" velocity. But the wave lengths are shorter in the 
rear and longer in front. If you follow a crest, say the 
last, you vill find that it is very quickly the next to last. 
AS now crestS appear at the rear, old crests traveling more 
slowly disappear at the front. This is an example of a single 
group speed but differing individual or “phase" speed. The 
phase speed, c, is related to length {for the pure deep water 


Sinusoid) by the expression 
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where T is the poriod and L is the wave length. 
Figure D-3 shows the observed frequency distribution 
ef apparent wave periods T (from Reference 16 in che form 


of a histeogran, 
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F = MEAN PEROD COMPUTED 
FROM SPECTRUM 
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Figure D-3 Observed Frequency Distribution of 

Apparent Periods of Waves in a Typical 
Irregular Sea 

If this distribution is measured at several times during 

a "steady" sea state, and each histogram is compiled from 

records of apparent period observed during a given period, 

and the increments T {in this case 1/2 sec.) are the same, 

the histograms will ail have approximately the same shape. 

Although this suygests the superiority of statistical over 

deterministic mathads, the frequency distribution does not 

provide enough information. 

Tt is spparent from Figure D-3 that observed frequencies 


are cantinuously diatributed Fram (in this case) about 


ase e ~y » 1.5. (T = 4 sees.) 
20 2 
“adh 2s YZ » 0,42 is ba 15 secs, } 


(The actual bandwidth is sumewhat larger, as wa shall see. 


The duts in Figure 2-3 correspends roughly to a gea state § 


or 6.) 
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A similar histogram could be derived from a single record 


aeons 
# 


such as Figure D-2, or a collection or ensemble of records, 
which furnish the distribution of wave heights. In the 
classical theory, however, there is no relation between 
wave period and wave height, except that the height-to- 
length ratio is limited to about 1/7 (Ref, 20). 

Joint distributions are equally unrewarding. What 


appears below is an example from Reference 22. 
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Feriad and Wave Hetght iu oe 42.7 ft/sec) 
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These histograms are for very small waves (2.9% max) that 
bos might appear on a pond or sheltered hay. Although the periods 
seem clustered about the average period, the lengths are very 


widely distributed about the “single wave" height. 
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A brief attempt has been made in this section to expose 

i the enormous complexity of sea waves. In terms of a possible 

{ mathematical descripticn, the sea is very irregular. Despite 
this, the conviction that a mathematical description is 
possible is greatly strengthened by the shsevvation that over 

a wide area, for periods of several huvrs, the sea may maintain 
an appearance which defies description in terms of "average" 


or “typical” wave lengths, or periods, but nonetheless appears 


to be steady. 


a Re “eb hen ee seth I ARS te. Agape An penta uettpebe Sanaa Pes Dla SRS, PA Sd Oni a 8 a is wie bot Me Beete cs Lae ne be oe AS Se Oo Ae eat Lan cE 


D-3 The Specification of the Irregular Sea 


4 y Numerous mathematical models for describing a real 


a 4 sea state have been proposed but all have proved to be 


inadequate except the description in statistical terms using 
the concept of the continuous energy spectrum, | 
It is very hepful in understanding the need for a 
Statistical description, to examine the deterministic models. 
It will be seen in each case that despite the deterministic 
character of the describing function in terms of specified 


amplitudes and frequencies, each method assumes information 


which can only be provided by statistical analysis of actual 
records. It should be noted that if any of the models 
described below were adequate, it could be used to provide 
wave forcing furctions in equation D~l, sinee the solutions ho 
for these equations are for a single periodic wave of specified 
amplitude and frequency. The complete solution for an aperiodic 
wave consists of the superposition ef an infinite for in the 
practical case, a very large) number of regulcr waves, in a 
manner which will be diacussed in the next section, 

Consider the wave recerd ef Figure D-$. 

(l) This may be represented by a periodic wave whose ' 
angiitude and frequen, has some significant relation to the 
irregular Wave. It is suggested by Sverdrup and Munk Ref. 22) 
that this can be done with a regular wave 
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wherein Xn and t are the average amplitude and period of the 


one third highest waves. The representation is shown in Figure 


D-5. 
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r(t} = tm cos(Ztt + ¢) = fq cosl(wt+ a) 
eit) w r(ter) 


Figure D-S Feriodie Wave Syaterm with an Amplitude 

Sompenent at a Single Spectral Frequency 

A serious sharteoning of this representation igs that about 

five sixths of the time the actual wave haights will be lower 

than the significant hoight of the idealized wave. A mere 

Sericus fault lies in che gelectica of the period +. A 

shige acta a8 a narrow band filter, responding oaly ta periods 

near the natuscsl period te If < 2&8 Boo much greater of less 

than i? Mery littie motion will reavlt. But the aetual aye 


record of figure G-$ contains wave components of 212 frequencies 
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and although the amplitude of a wave of period +t ® T, may 


be quite small in comparison to r_, its effect on the motio. 


m!’ 
of the ship may be quite large. The idealization of "i yure 
D-5 has filtered out these components. 

(2) The irregular wave could be represented by a Fourier 


series with various amplitude components at many discrete 


spectral frequencies. This representation is shown in 


Figure D-6. 
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Figures B-& Periodic Wave Syseten with Amplitude Cemporents 
at Many Uareratée Spectral Frequencies 
ft if based Gn a harmonic analysis of the actual wave revord 


Of Figure Geb Tor the gerict t,., and For that pexiod only is 
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be made as exact as desired by taking N as large as necessary. 
But there is a severe practical limitation. In actual wave 
records, the amplitudes of harmonic components are neglible 
for periods less than 25 seconds which means that for n < (t,/25), 
all terns would drop cut. The harmonics of maximum amplitude 
will have a period which is greater than the waves of siqnificant 
(one third highest waves). The amplitudes of higher harmonics 
with periods less than 3 seconds would be neglible. Thus 
LE t, were chosen as 
t) = 400 szconds 30: minutes 
Neglected lower 
harmonics, n= 1, to 4 70 
Neglected higher 


hermorics; n > 30 §00 


The significant portion of the Fourier series for the 400 
second case consistsof about 36 ~ 4 3 26 texms, and fer the 
30 jsinute cage of about 600 ~ 70 » §30 terme. Now it ia 
@lse true that lew eetual waves in the Fourier series repre- 
gantation age cauged by the phase cancellation ef large 
nucbers ef ba: nie componente of shal amplitude and laree 
actual weaves are Caused By phase reinforecsent af the same 
Cetponenis. The phase angles é, ¢an be determined exsesly 
by harmonic analyais for the given pericd t 7” pe hut will 


agsecssarily differ for ony other tise posed. 
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(The temptation for engineers to think in terms of this 
modi] is unbelievably strong. After reading to the end of 
this appendix, the reader should consult Reference 23. In 


this paper, the author accepts the inevitabilit. of a statisti- 


cal description ‘n terms of the energy spectrum, but 


reverts to a function of the form shown in Figure D-5 to 
approximate the spectral @istribution!) 
(3) A third possibility for the description of actual 


wave recc~ds is by use of the Fourier integral. Such a repre- 


a - sentation is shown in Figure D~7. 
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‘ £"; In this manner r(t) can be made to represent the actual 
= ree : 
: wave record over any length of time chosen for analysis and 
B 
: this representation would be exact for the interval analyzed. 
: It would fail, however, for time outside the interval in which 
e 
t the analysis was performed, since the record would either have 
A to be defined as identically equal to zero or would remain 
: unknown. 
One disadvantage of this method is that there is no known 
; precise procedure by which, starting with a wave record, it 
: is possible to determine the appropriate a, (u) and b, (a) 
Nevertheless, such a procedure is theoretically possible. 
The Fourier Integral method would be a convenient one for 
studying proble sn transient response. It does not, 
however, lend itsclf to problems of the seaway because of 
the extrumely long duration of the records to be analyzed. 
The three modole discvssed above replace an actually 
measured wave profile of sore definite leagth with a harmonic 
i approximation. The most serious fault in gach method is that 
the informution ‘termined in each case applies only ta a 
q Specifiod tima interval. The essential irregularity of 
4 thu saa ia never adequately acceunved for, although the 
4 aeproxination to a given racord could theoretically be made 
d ag ¢loae as desired. In additien thera ares in wach case 
4 prectical difficultics which have boen deseribed whie 
= | sivlously vitiasc the usefuliness of chase methods. 
1 4 However, for any given moa state, over appreciable 
\ : lengths of tine, as santioncd abave, thera is a “steady 
q : | D-16 
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state" appearance of the sea. This has led cvceanographers 
to make use of che concept of the energy spectrum, exploiting 
the statistical representation of a random process in a form 
similar to that developed in communication theory by Tybey and 
Hamming, Reference 17 and Wiener, Reference 18. 

The motive in discussing the three mavhematical models 
described above is not only to exhibit their inadequacies 
and the need for a random statistical approach as will be 
described below but to clarify a loss of detail which must 
be accepted as a consequence of the use of the statistical 
method. Although in every way favorable to the “deterministic" 
models described above, the energy spectral method requires 


that the phase angles are random and independent. 
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D-4 The Statistical Representation of the Random Sea 


According to this representation, the wave amplitude 


‘ata point can be represented by the integral 


r(t) = ie cos [wt + e(u))¥ Is. (w) aw D=5 


+t will not be necessary to use this integral (which is not 
an integral in the oréinary sense). This integral, which 
can be approxinated as closely as required by a sum of dif- 
ference equations, is distinguished from the previous de~ 
scribed models by the random nature of the phase shifts e(,)) 
which are equally prebable and equally independent. 

The most important term in equation D-5 is the spectral 
eneruy Sw). The meaning of this term may be visualized 


by examining Figure D-8. 
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Here it is assumed that the typical wave record r(t) is 


composed of a large number of sinusoidal components and 


that the alevation of the sea surface at any time t may be 
. a represented by the sum 
r(t) = lim 2 cos [w, + e(w)) [25, (w) 60) 1/2 D-6 
Or 
GW? oO 
The term S,(w), a continuous function of wave frequency, 
and is referred to as the spectral density. It comprises 
the totality of the statistical description of the sea 
state. In Figure D-&, it can be seen that S.(w) is a 
function which fer any increment $y is proportional to the 
total wave energy in that increment. In fact for the 
increment Su at the central frequency We the total energy | 
is 


Pg [S. (u) Su} a-7 


where p is the densit, and g is the acceleration constant. 
Tha total energy of the wave systam rcepresanted by all the 
component! energies in the spectrum is given by the integral 


rg f & (.}du 
rg b ! ree! 


This integral, which reprosgenta the area under the spectral 


curve of Figure Ge@, is unually called E and hence the total 


y+ 
é 


energy iB ogf.S, haw the units of (Leagth” x time). 


Wihie the energy apcctrum $b) has the appearance af 


an arpiirudy socetrun, it du different in a very significant 


Way Vick Oxplaing Lts unique ability to characterize the randan 
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sea. The potential energy of a single sinusoidal wave (per 


unit length of crest) can be represented by 


1/2egr? 


where r is the wave amplitude. It follows, then that at 


a given central frequency w a fictitious wave having the 


n! 
same eneryy as represented by equation D-7 would have the 
height 

1/2 


ns [2S. (w)de] 


The spectral density S,(w), however, represents the cumula- 
tive effect of all measurable amplitudes in the increment 
Sw centered at Wes 

S,(s) can be readily determined from amplitude time 
histories such as Figure D-3. The energy spectrum is in 
fact the Fourier transform of the autocerrelation function 
of a given amplitude time history. Thus 


r(t) réter) o*ST ar D-8 


’ L 
S, (0) = oy f 
where t iS a continuously varying time lag. 

There are two cardinal paints to keep in mind in asses- 
Sing the value of the spectral density methed. The funetion 


~ 


g(t) given by equatieas De5 and D-G, although it represents 
&@ statistically Valid profile will net be used in the subac- 
quent develognent al shin metions. Gnily the apeceral density 


Bots) du required, aAs will bo shown in the nexe aectien, 


thas appreuach pravites «. stulistieal deseription oF ship 
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Figure D-3 ‘Yypical Sea Spectra 


motions in the frequency domain. ‘The implication of this, 
Whicn has been alluded toa abave, is that this approach 
discards all information which relates to the relative 
phase of ths varicus frequency components. The affect 
af this will ba to introduce a deyree of indaterminacy 
in the equetions af gireratt mokian whose :mpertance wiil 


be vory hard ta assess. 
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D-5 Ship Motions in an Irregular Sea 


The solutions of equation D-1 for each frequency 
are linear solutions of the form D-2. I1f we knew the 
harmonic components of each frequency band é in Figure 
D-9, specifically the amplitude of each component Ws 
then we could combine the motions by the principal of 
superposition and determine a total motion which applies 
to all forcing functions present in the sea pattern, 
This, however, requires that we must use one of the 
three deterministic models shown in Figures D-5, D-6, 
and D-7. These are, in fact, the deterministic medels 
which cannot be used for the reasons discussed in the 
previous section, The use of the sea energy spectrum 
ef one of the experimental forms shown in Figure D-9 
Fequires that we use the statistical information which 
these spectra cantain in a tnmanner which is within the 
applicability oF the principal of sumerposition, 

The stetistical sea ensrgy spectral data is used 
with the solutions ef equation D-b in the following way: 
A spoctraum of ship response amslitudes es a fenecion of 
(encounters) frequency is dotertmined by the use af the 
enetdy spectrun Sy Ge) and the ship response: Gporater 
ele) by the Formula 

Sy Coed a plu} Sy t+) b-$ 
The ship respense oper.tion is defined as the squared valus 


of the respensr tp GAAS) (sce equation D-2) 
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r(w) = [ns (w)/Ay (w))* 
3 J 


for each significant value of w. Using the symbol 
([R(w)]? for r(w), equation D-9 may be written 


S_(w) = [Ru]? + Sy (w) 


a fori in which it often appears in the literature. (See, 
for example, p. 31 of Reference 2, and page II-8, gic eeneee 
The individual results in the calculation of equation 

D-9 are illustrated in Figure D-10 and will be explained 
in some detail below. Figure D-10 shows the application 
of the principle of superposition. for sea spectrum at a 
fixed point (Figure D-10(a)) for two typical ships: a 250' 
ship moving at 7.97 knots (Froude number Fr = 0,15) and a 
§00' ship moving at 11,27 knots @r = 0.15}. In terms 

of the frequency of encountar, », the eneryy spectra are 
shifted ta the right by the transformation 


ul ty 
oy = {l a oo 


where wy 38 the wave frequency and U is the ship speed, 
In addition, a@ explained in Refersnee 2, the spectral 
amplitudes are transfarmad by multiplying each amplitude 


by the Jacobean of the transformation: 


i 


ie S24 

g 
Naturally, tha total areas under the original and transformed 
curves ara identical. In Figure 0-108 the curves are smoothed 
curves drawn through a large collection of individual 


Balutions to cquatinn D+il. For example, at the center 
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frequency ou, = 0.8, Figure D-10b shows that the ship 
response in pitch, scaled by the wave amplitude A. is 


about (for the 250' ship), 


(8/5)? =2x i074 (= Ng/Ag) ” (Continued page D-25) 


meet 


Notes for Figure D-10: 


(1} The units of each of the sets of curves are as 
follows: 


Fig. 9(a) 2 


i 


Units of S$. (w) are FT“ ~ SEC. 


1 


it 


Units of 6/5 are (rr? - SEC?) ~ 
1 


Fig. 9(b) 


Fig. 9(c) 


i} 


Units of S,(w) are (SEC) 


(2) The spectrum &.(w) may either be determined 
numerically from an actual wave (one-dimensional) record 
such as that in Figure D-2 or anelytically by a formula 
such as the Pierson-Moskowltz equation cited on p. 25, 
Reference 2, ar the Neuman Spectra (equation 1.28, Reference 7). 
In the fermer case, the auto correlation of the wave record 


Ls obtained 


o(v) = him fo ale) a(t + 1) dt 


The Fourler trangterm of 2300) ig the power spectrum for 


the Variable z, 


st eugene are} 
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In the lattar case, the amplitude ef the sea Spectrum 
is expressed For any given Traguency 28 & funetion ef sone 


parametor, for onample @igniisieant wava haight bi {in the 


Plarson Moskowitz cguation). 
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Notes for Figure D-10 - Continued 


2 oe -2, 2 
, where A= (8 x 10") g ' 
B= , h 
(33.56) 1/3 
w = wave frequency 
or (for the Neuman spectrum) as a function of wind velocity 
bh 
. 2 
S,, (w) = (c/w®) exp (-297/t 7?) cos’ xm 
; c = 32.9 FY - SEC? 
U = wind velocity 
i angle between wind and wave patterns. 
ket “4 
; wore (} RESPONSE. 
a6) Sena seca 7 ‘eters 
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3 where € (and ",) is the pitch in radians and A (and nj) are 
3 i 
As 
the wave amplitudes, and about 1/2 that fer the 580' aiip. 
4 
: 
The spectrum el Figure h-lOb le the ship response amelicude 
4 operator ruferrud to in equatien O-98 above.aa r(e). (Note: 
q in the literatur.s the various ship response amplitude opere- 
4 ; 
4 terse ate defined and nen-dimensionalized in various ways. 
: For exatpla, in Reference 2&2, the response amplitude opera- 
ter (RAD) in pitch is scaled by 2x wave amplitude/ship length: 
i D-25 
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@L/2S, L is the ships length and g is the wave smplitude, 
In the case of that reference, equation 5 is used. The 
difference in the curves for the 250' and 500' ships are 
very significant. They illustrate the general principals 
that apart from hull shape and weight distribution, that 


wave components of only certain frequencies have significant 


effects on ship motions. In general: 


(a) Wave components cf about 3/4 ship length or above 
have appreciable effects in pitching motions in head seas, 


{o) The most severe motions result at about the ships 
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natural pitchi:g frequencies, The same kind of statements 


can be made about all Ship motions components. i 
The tetal pitch ship response in a particular sea state ' 
i 

represented by the sea state spectrum of Figure 0-104 is 


aro 


er determinad by multiplying ordinate by erdinate the apectra ; 
' 4 of Figure O-]0a and D-lGb ta obtain the smoctrum of Figura 
s q D-lde. Adain, it will he noted that there is a significant 
q q chanve in the shape ef the spectra. For exanple, the pesk 
7 4 amphituda of the RAO far the 240‘ ship is at about s fra- 
q . q qucney ef 1.2, Khareaa the piteh ampiituds spectra fer the : 
2 q Same shie is velativaly flat between frequencias ef 0.6 
q 4 aod lie radians/second, The meaning of ihe shifting peak 
: 3 amplitudes in the three diagram’ should be clearly vuaderstoom. ' 
s ’ . ia) In the sea apoctrom Fievre G-lOe the peak arpld- } 
: tude of the original spetctrus is at about o-* 9.4. The ; 
forward apecds of the two chips shifis the Frequency of en- 7 
counter slightly to the vight since the ships sca directed | 


a wales as ase 
en ont. 


into head-on seas. 

* (b) The smooth curves of ‘igure D-i0b represents an 
interpolation through an ensemble of solutions to the pitch 
equation of equation D-1, wherein the pitch @ has been 
scaled by the wave amplitude. That is 

(0/0)? = (ne/al)? 
where ao is the amplitude of the exciting wave. Each 
point on the smoothed RAO curve of Figure D-l0b represents 
a Single independent solution of equation D-1 (in this case 
the pitch equation). That is to say, the point on the 
raired curve (ut)? quoted above at an encounter frequency 
of 0.8 (for the 250' ship) is the Ng (ex pitch) selution to 
equation D-] for the (unit) wave foreing funetion of 


z/a, = cos (0.8%) {Read-on) 


Likewise, the maximunm amplitude fer the 250' ship which 
OCGUrK at an encounter frequency of about w » 1.1 radiana/ 
sesond is the eingle particular ars solution for the particular 
{unit} wave 
afa., mcoes 2.168) 

For cach metian equation Bel most be ealved to deterring 
Separate zy each point on the Rav. In Reference 2, For 
exstpie, Solutions fer pateh and heave sre calculated at 
increments oF the encounter irccuensy of 0.2 from 1.0 to 6.4. 

{o} The pitch amplitude spoctra Figure O-l0e is fer 
each freguency, the proluet of the transformed sea spectra 
fin tease of encosnter Frequcncy) and the RAO of Fisure beldb, 


However, what 49 important to realice is that the piteh 
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Spectrum (89 (a) of Figure D-l0c is itself a spectrum of 
exactly the same statistical nature as Figure D-10a. The 
Spectral curve tells us nothing directly about the expected 
amplitudes of pitching and nothing whatsoever concerning the 


frequencies at which these expected amplitudes might occur. 
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i Statistically, the expected amplitudes are related to the 
| total area under the pitch amplitude spectra, Figure D-luc. 
Thus for tne total energy 
! Ey = f Sytw)du = 1/2 fe 2 Sg} du 
4 statistical thcory states that for a typical “Rayleigh” 
: ; distribution, the average amplitude Bae “significant” (or 
q average of lf: largest ) Jy 731 and average of rhe 1srot® 
q , . highest. amplitudes sro! are given by the expressiuns 
q : os 8 1.253 EL 


re) 
5) 
210” 2.546 Ee 
Thus, in the CREE qgucéced in Figure Pel, the tetal “energy” 
in tha curve for the 50' ghip From aswe © is about 


°% = 
agente? tweeltee@ic! anus, 


x, 
% 
+ 


+ 3,8" 


ad 


4 a ; . 
2a) 2.9 


Uy cee Sas 8" 
iia = 
Althuugh this fect has beun veforred to with varying tegrecs 


oY erphasia @averai tinea in the Eeregeing discussien it 


ls comegavlersd viryvy impnrta.t in View oF what witli fol low 


to pnlaras upon it agein: 


(where © 


Although the various "statistical" amplitudes referred | | 
to above ("average", "significant" or "1/10°" highest") ave Gy 
functions of the total area of the energy Curve, this method 
provides no information to associate these amplitudes with . 
any specific frequency. Although our mathematical intuition | 


assures us that there is no other way for it to be, 


it is 
nonetheless a feeling of somehow beiny cheated that everosmes 
us when we review the three sets of curves in Figure D-10. 


The curve of Figure D-10b /s composed of a faireg interpo- 


lation of points, each of which as a solution of one of the __ 
equations of the set equation D-1 is a solution which com- 


bines both amplitude and frequency (as well 38 phase lag) 


information. For example, aS quoted above, at an encounter” 
frequency of 0.8, where (8/1) m2 xX 16°", this Ship pitches 
about - a 

0 = 0.815°/FT of wave | 
That is to say that for a 2' wave of the fora 


a2. ¢ tos 0.et 


The ship would pitch according to the axpresaian 


wn 
i 


or 


sy 4 


.* 2.815 (eos G.8t * £5) 


FI although determined by the solution vf equation: D-2 


is not included in the data of Figure D-10D). However, in 


the gineration of Figure D-l0e, both “he frequueney and phase 
information of solution De? is discarded, and amplitudes of 
given Ginengions (such as “average” or “signifacane® o¢ 
“average of right highost") or the probebility of an 
ampiltudia of any givun height Jo dutermined from the cota 
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integrated area of the S,(w) curve. For example, statistical 
theory for this type of distribution (the "Rayleigh" distri~ 
bution) states that the probability that the amplitude 8, a 
given response will exceed an amplitude €, given height Ore 
is 

P(0>6,) = exp (-02/2m) 
Since the “significant” amplitude “average of the 1/3rd highest) 
for a Rayleigh solution is 


8 = 2.00 E 


1/3 ) 


‘where 


E = ; Sy (w) qs 


fox any sufficiently large number N such that N = P (vy? 4g) 
the ost probable larges: pitching amplitude 2x is 

ay zs By ya Inn / 2 
No similar information is provided from Figure D-1l0c. for tne 
most probable fraquency or range of frequeacies for which a 


given amplituae may occur. 
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APPENDIX E 


COMPILATION OF HARRIER DATA 


COMPILATION OF HARRIER DATA 


The data on the Harrier GR Mark 1 aircraft which follows was 
received by CADCOM after the conclusion of the project described 
in this report. It was furnished by Dr. Bernard H. Carscn of the 
Office of Naval Research, London. 

CADCOM includes it herein eerietiy for information purposes 


and makes no representations about its accuracy or completeness. 
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